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We present a quantitative analysis of the phonon lineshapes obtained by neutron inelastic scat- 
tering methods in the relaxor ferroelectric Pb(Mgi/3Nb2/3)03 (PMN). Differences in the shapes 
and apparent positions of the transverse acoustic (TA) and transverse optic (TO) phonon peaks 
measured in the (300) and (200) Brillouin zones at 690 K are well described by a simple model 
that couples the TA and soft TO modes in which the primary parameter is the wave vector and 
temperature-dependent TO linewidth T{q,T). This mode-coupling picture provides a natural ex- 
planation for the uniform displacements of the polar nanoregions, discovered by Ifirota et al, as the 
PNR result from the condensation of a soft TO mode that also contains a large acoustic component. 



I. INTRODUCTION 

Considerable advances have been made over the past 
several years in our understanding of the lattice dynamics 
of the relaxor ferroelectric Pb(Mgi/3Nb2/3)03 (PMN), 
which is a prototypical member of a class of compounds 
that possesses exceptional piezoelectric properties, and 
an unusually broad and frequency-dependent dielectric 
susceptibility that Maks at a temperature Tmax (~ 265 K 
at 1 kHz for PMN) B Pioneering neutron scattering work 
on PMN in 1999 by Naberezhnov et aln has been fol- 
lowed by a systematic seriea_|Q£_neutron scattering 
periments by Gehring et alBuWB and Wakimoto et alU 
The results of these latter measurements have firmly es- 
tablished the existence of significant soft mode dynamics 
in PMN. In particular, they have demonstrated a linear 
temperature dependence of the zone-center soft phonon 
energy squared (hujo)'^ over a wide temperature range, as 
shown in Fig. 1. 

A fundamental feature that appears common to all re- 
laxors is that of the polar nanoregions, or PNR. The 
presence of these regions was inferred in 1983 by Burns 
and Dacol, who found that the optic index of refraction 
of PMN deviates from a linear temperature dependence, 
at a temperature Td « 620 K that is far above Tmax-u 
In conventional ferroelectrics this deviation corresponds 
to a uniform polarization that only develops below Tc- 
Given the absence of any net polarization in PMN at Td, 
the optical data indicate the formation of tiny regions of 
local and randomly-oriented polarization that were spec- 
ulated to originate in Nb-rich parts of the crystal several 
unit cells in size. The dynamical effects of these PNR 
are clearly manifest in neutron scattering measurements. 
Gehring et al. have found that the long- wavelength trans- 
verse optic (TO) phonons are heavily damped below 
which gives rise to the anomalous "waterfall" feature 
Neutron scattering measurements further suggest that 
the appearance of the PNR at Td is accompanied by a 
remarkable overdamping of the zone center (q — 0) TO 
phonons. As shown in Fig. 1(a), the zone center TO 



phonon linewidth broadens with decreasing temperature 
until finally at Td no discernible phonon peak remains. 
It has subsequently been reported that the ovecdamped 
zone-center TO mode reappears below ^ 220 K. El 

The formation of the PNR at Td is the most impor- 
tant aspect of the relaxor problem, and it is intriguing 
that it is accompanied not only by the overdamping of 
the zone-center TO soft mode, but also by the appear- 
ance of diffuse scattering. The Brillouin zone dependence 
of the neutron diffuse scattering intensitiea-iias been ex- 
tensively measured by Vakhrushev et a/., □113 and later 
it was shown by Naberezhnov et al. that the onset of 
the diffuse scattering occurs at TdB The coincidence of 
these phenomena at Td point to a picture in which the 
PNR result from the condensation of the soft TO mode. 
If this picture is correct, then the atomic displacements 
associated the diffuse scattering must match those associ- 
ated with the soft TO phonon, namely they must satisfy 
the center-of-mass condition. However, the diffuse scat- 
tering intensities in different Brillouin zones are entirely 
inconsistent with the TO phonon intensities. This longs 
standing puzzle was resolved recently by Hirota et al.t3 
By calculating the dynamic structure factor with the as- 
sumption that the soft TO mode condenses with an addi- 
tional uniform shift of all atoms, they demonstrated that 
the intensities of TO phonon and diffuse scattering can 
be completely reconciled. Their calculation shows that 
this uniform shift 5 shift is comparable to the displace- 
ment due to the TO mode condensation Sc.m., and that 
S shift is 60 % of the displacement of Pb atom. While 
the origin of this uniform phase shift is unknown, the 
success of the model provides further support to the idea 
connecting the PNR and the soft mode. 

In this paper we demonstrate that a simple model that 
couples the transverse acoustic and soft optic modes, sim- 
ilar to that used by Harada et al. to describe, the asym- 
metric phonon lineshapes found in BaTiOaJld is able to 
reproduce the phonon lineshapes observed in PMN in 
both the (200) and (300) Brillouin zones, which exhibit 
very different scattering profiles. This naturally implies 
that a condensation at Td of a coupled TO mode would 
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FIG. 2: Temperature dependence of the (2, 0, 2) Bragg peak 
linewidth measured along the (a) [101] (longitudinal), and (b) 
[101] (transverse) directions. If the system transforms from a 
cubic to a rhombohedral phase, then the (2, 0, 2) peak should 
split along the longitudinal direction. Instead, only the trans- 
verse scan shows a clear change near Tc. For ease of compar- 
ison, panels (c) and (d) show the (2, 0, 2) lineshape measured 
along the longitudinal direction at 300 K and 50 K. The solid 
bars represent the instrumental longitudinal g-resolution. 



FIG. 1: (a)Temperature dependence of the (2,0,0) zone cen- 
ter TO phonon linewidth Fi . Above Td Fi increases gradually 
with decreasing temperature, whereas below Tc it decreases 
rapidly. For temperatures between Tc and Td it is impossible 
to determine the linewidth due to the overdamped nature of 
the TO mode. Solid lines are guides to the eye. (b) Tem- 
perature dependence of the (2, 0, 0) zone center TO phonon 
energy squared (hujo)^ from Ref. 7. The dotted line shows the 
Curie- Weiss behavior of the dielectric constant 1/e cx T — To 
with To = 400 K reported by Viehland et al. 



then contain a TA component, and this might closely 
reLate to the uniform phase shift reported by Hirota et 



II. EXPERIMENTAL DETAILS 

Single crystals of PMN were grown using a top-seeded 
solution method with PbO as flux. Two crystals were 
chosen for the neutron scattering experiments. One has 
a volume of 0.09 cm^ and a mass of 0.74 g, and the other 
has a volume of 0.40 cm'^ and a mass of 3.25 g (the density 
of PMN is 8.13 gm/cm'^). Since both samples show quan- 
titatively similar results, we will not distinguish between 
the two samples in this paper. 

Our neutron scattering experiments were performed at 
the BT9 triple-axis spectrometer located at the NIST 
Center for Neutron Research. The (002) reflection of 



highly-oriented pyrolytic graphite (HOPG) crystals was 
used to monochromate the incident neutron energy Ei 
and to analyze final neutron energy Ef. The data 
were taken with a fixed final energy Ef — 14.7 meV 
(A = 2.36 A^), and horizontal coUimations of 40'-40'-S- 
40'-80' ("S"=sample) between source and detector. An 
HOPG filter was placed after the analyzer to eliminate 
higher order wavelength contamination in the scattered 
beam. The crystal was placed with a natural {100} cubic 
facet facing down on to a boron nitride post, and held 
in place with tantalum wire. This orientation allows ac- 
cess to (/i0Z)-type reflections in the horizontal scatter- 
ing plane. The mosaic spreads of the crystals were less 
than 24' at the (200) reflection which is resolution-limited 
with the above configuration, and indicate a high crys- 
tal quality. The room temperature lattice constant is 
4.04 A, therefore one reciprocal lattice unit (rlu) equals 
1.553A-1. 

PMN reportedly remains cubic down to 5K.Ej Con- 
sistent with this observation, we find no definitive evi- 
dence of a structural transition in the temperature range 
50 < T < 400 K. On the other hand, a macroscopic fer- 
roelectric state can be established in PMN by cooling 
it in a moderate {E — 1.7kV/cm) applied electric field. 
The induced polarization then vanishes upon warming 
above Tc = 213 K as a first-order phase transitionjj To 
date, the aforementioned reappearance of the zone-center 
TO mode, which also coincides with the disappearance of 
the TA phonon broadening ,0 is the most direct evidence 
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of ferroelectric ordering in zero field associated with Tc- 
Other anomalies have also been reported near T^, such as 
the sharp peak in the temperature dependence of the hy- 
personic damping reported by Tu et al. ,tj and the abrupt 
change in the thermal cpqpansion of the lattice parameter 
observed by Dkhil et aZJlj To pursue this issue more care- 
fully, we looked for subtle structural changes near Tc by 
employing a tighter sequence of horizontal coUimations 
(40'-10'-S-10'-80') to improve the instrumental q resolu- 
tion. Figure 2 shows the temperature dependence of the 
(2, 0, 2) Bragg peak width along the (a) [101] (longitudi- 
nal) and (b) [101] (transverse) directions. If PMN were to 
transform from a cubic to a rhombohedral structure, then 
the (2, 0, 2) Bragg peak would split along the longitudinal 
direction. However, we observed neither a splitting nor 
a broadening of the (2, 0, 2) Bragg peak along [101]. Fig- 
ures 2(c) and (d) show longitudinal scans of the (2, 0, 2) 
Bragg peak measured at 300 K and 50 K. The peak width 
is almost resolution limited for 50 < T < BOOK. Instead, 
we observed a small, but clear, jump in the transverse 
width at Tc which is shown in Fig. 2(b). However, we did 
not observe a corresponding enhancement of the peak in- 
tensity below Tc that is normally expected at a structural 
phase transition due to the release of primary extinction. 

III. MODE COUPLING MODEL CROSS 
SECTION 

The neutron inelastic scattering technique provides a 
direct measure of the scattering function ^(q, w), as this 
is simply related to the experimentally measured scatter- 
ing cross section. More importantly, 5(q, lo) is related to 
the imaginary part of the dynamic susceptibility x" (q, ^) 
via the fluctuation-dissipation theorem 

5(q,c.)-{l + n(c.)}x"(q,c.), (1) 

where n(uj) is the Bose factor {^e^/f^sT _ l)^i. By choos- 
ing the scattering vector Q = G + q, where G is a re- 
ciprocal lattice vector and q is a phonon wave vector, 
one can measure the phonon dispersions in different Bril- 
louin zones. The measurements reported here were made 
at G = (2,0,0) and (3,0,0). However a model cross 
section is required in order to extract meaningful param- 
eters from the measurements, and to correct for the ef- 
fects of the instrumental resolution. The_£orm we chose 
for 5'(q, w) is discussed by Harada et aZ.,E3 and pertains 
to a system in which a coupling exists between two vi- 
brational modes. In this case we have 

5(q,a;) = {1 + n{uj)] {[{^1 - u;^)B - T^AlFi 
+2XBF1F2 + ml - u;^)B - riA]F^}, 

(2) 

where 

A^{nl-u^){nj^uo^)^u^TiT2~X^, (3) 

B = Ti{nl-uj^) + r2{nl-uj'^). (4) 



The model parameters ft, T, F, and A represent the 
phonon frequency, damping, dynamic structure factor, 
and coupling constant, while the indices 1 and 2 denote 
the TO and TA modes, respectively. Hereafter, we shall 
refer to this as the mode-coupling (MC) function. This 
function shows larger coupling effects the closer the TO 
and TA frequencies become, or the more heavily damped 
either mode becomes. In the appendix we identify several 
important features of the MC function that can produce 
large differences in the phonon profiles. 



IV. NEUTRON SCATTERING 
EXPERIMENTAL RESULTS 

A. Phonon profiles at 690 K 

Substantial differences between the phonon profiles ob- 
tained in the {2,0,q) and {3,0, q) Brillouin zones are 
apparent in Fig. 3, which displays typical data sets for 
q ~ —0.2 measured at 690 K. A major difference is that 
the scattering associated with the TA phonon appears to 
peak at different energies, i. e. the TA energy seems lower 
at (3, 0, -0.2) than it is at (2, 0, -0.2). This difference is 
systematically observed at different q positions as shown 
in the inset of Fig. 3, where the apparent TA phonon 
energies obtained from the peak position are plotted as a 
function of q in the (200) and (300) zones at 690 K. The 
solid lines are merely guides to the eye. As is clearly seen, 
the TA mode energies in the (300) zone are consistently 
lower than those in (200) at all values of q studied. Of 
course, for a crystal that possesses true long-range trans- 
lational order, all phonon energies in different Brillouin 
zones must be identical. 

When fitting these phonon profiles with the MC func- 
tion, we found that a change in the sign of the product 
of the dynamic structure factors F1F2 can reconcile the 
discrepancies in the TA phonon peak position. Specifi- 
cally, if we set F1F2 < for (2,0,-0.2) and F1F2 > 
for (3,0,-0.2), we obtain the same TA phonon energy 
5I2 for each zone. This occurs because for positive F1F2 
the spectral weight of the TO mode increasingly shifts to 
the region below the TA mode energy 1^2 with increasing 
TO mode damping Fi, whereas for negative F1F2 this 
spectral weight shifts to the region between the TO and 
TA mode energies. Thus the MC function produces a 
higher (lower) apparent TA phonon energy for negative 
(positive) F1F2, even though ^2 is always the same. Fur- 
ther details are given in the Appendix. The fits to the 
MC function include an additional Gaussian component 
which is used to describe the elastic incoherent peak at 
huj = 0. The fitted curves are shown in Fig. 3, where the 
solid and dashed lines correspond to the MC and Gaus- 
sian functions, respectively. The fitting parameters are 
listed in Table 1. It is satisfying to note that the values 
for and F derived from the fits of both the (2, 0, —0.2) 
and (3, 0, —0.2) phonon profiles agree quite well, even 
though the apparent TA peak positions are different. 



4 



PMN : T = 690 K Ef=i4.7mev, 40-40-40-80 



:d 
o 
o 



CO 




4 8 12 

■Rco (meV) 

FIG. 3: Constant-Q scans measured at (a) Q — (2, 0, —0.2) 
and (b) Q = (3, 0, -0.2) at 690 K. Solid lines are fits to the 
mode-coupling function convoluted with the instrumental res- 
olution function. A negative F1F2 is used for (2, 0, —0.2) while 
a positive F1F2 is used for (3, 0, —0.2). Although the scatter- 
ing intensities from the TA modes peak at different frequen- 
cies, the fits yield the same value for ^2 ~ 4.2 (as is required), 
as well as the same value for A = 17. Note that the TA struc- 
ture factor at (3, 0, —0.2) is « 100 times smaller than it 
is at (2,0,-0.2). The inset shows the apparent phonon fre- 
quencies of the TA modes determined solely from the peak 
position. The TA modes at (3, 0, q) appear systematically at 
lower frequency than do those at (2, 0, q). 
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FIG. 4: Constant-Q scans measured at (a) Q = (3,0, —0.12) 
and (b) Q — (3,0,-0.08). Solid lines are fits to the mode- 
coupling function convoluted with the instrumental resolution 
function. Here F1F2 was chosen to be positive for both scans. 
All parameter values are listed in Table 1. 



the actual intensity of the TA mode at (3, 0, —0.2) (which 
is proportional to |F2|^) should be over 100 times smaller 
than that at (2, 0, —0.2). Yet the TA peak is clearly visi- 
ble at (3, 0, —0.2). This can also be explained by the MC 
function, since the intensity of the TO mode is trans- 
ferred to the TA mode due to the coupling, which in turn 
produces an enhancement of the TA cross section. This 
feature is also described in more detail in the Appendix. 



Another important difference between the (200) and 
(300) zones is that the dynamic structure factor IF2I 
at (3,0,-0.2) is about ten times smaller than that for 
(2, 0, —0.2). This fact is readily confirmed since, for wave 
vectors close to the zone center, the ratio of jFaj in two 
different zones should agree with that for \Fn\, the nu- 
clear static structure factor. This is given by 

^^Jv(Q)-5]6,e^'^•'^^e-^^ (5) 

i 

where bj, dj, and e~^^ are the neutron scattering length, 
atomic coordination vector, and the Debye- Waller factor 
for the jth atom, respectively. In fact, this formula gives 
|FAr(200)|/|FAr(300)| = 11.2 which is in excellent agree- 
ment with the value |F2(2, 0, -0.2)|/|i^2(3, 0, -0.2)| = 
12.0 obtained from the fits. This implies, however, that 



Features similar to those discussed above are also found 
in the phonon profiles at smaller q. Data taken at 
(3, 0, -0.12) and (3, 0, -0.08) are shown in Fig. 4(a) and 
(b), along with the fits to the MC function. As before, 
the solid and dashed lines correspond to MC and Gaus- 
sian functions, respectively. The fitting parameters are 
listed in Table 1. The coupling constant A was deter- 
mined such that r^i and 1^2 are equal in each zone. At 
smaller q, the parameter Fi becomes larger, meaning that 
the TO mode is more heavily damped. This is consistent 
with the "waterfall" picture. We have thus shown that 
the phonon profiles in both the (200) and (300) Brillouin 
zones can be explained satisfactorily by the MC function 
at all q. However there is still an open question. While 
we know that A decreases at smaller q, we do not yet 
have a specific model that describes the q-dependence of 
A. 
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TABLE L Parameters obtained from the fit of the 690 K data 
to the MC function. 



Q 




Qi 


^2 


Ti 


Fa 




1^2 1 


A 


(2,0, 
(3,0, 


-0.2) 
-0.2) 


10.85 
10.92 


4.30 
4.20 


6.10 
6.05 


0.79 
1.95 


27.0 
27.0 


13.1 
1.09 


17.4 
17.4 


(3,0, 
(3,0, 


-0.12) 
-0.08) 


9.05 
8.76 


2.67 
1.92 


7.77 
13.1 


0.69 
0.03 


25.7 
31.3 


< 0.5 
0.64 


7.21 
1.88 



TABLE IL Parameters obtained from the fit of the data in 
Fig. 5 to the MC function. 



Q 




T{K) 


Qi 


^2 


Li 


La 


\Fi\ 


1^2 


A 


(2,0, 


-0.2) 


700 


11.1 


4.60 


5.24 


1.24 


10.5 


5.65 


17.4 






500 


11.1 


4.45 


6.07 


1.70 


12.5 


6.07 


16.8 






290 


12.0 


4.90 


5.96 


1.52 


11.5 


6.04 


17.9 


(3,0, 


-0.2) 


570 
350 


11.2 
11.2 


3.90 
4.29 


6.00 
6.00 


2.89 
4.81 


12.1 
14.0 


0.78 
0.82 


17.4 
17.4 



Temperature dependence of the phonon cross 
section 



To check the consistency of our model, we have an- 
alyzed the TA phonon profiles at (3, 0, —0.2) at dif- 
ferent temperatures. Before showing the results for 
(3, 0, —0.2), we first present the temperature dependence 
of the phonon profiles obtained at (2,0,-0.2). These 
are displayed in the left panel of Fig. 5 at (a) TOOK, 
(b) 500 K, and (c) 290 K. Each profile has been fit to a 
resolution-convoluted MC function plus a Gaussian cen- 
tered at ?iw = with F1F2 < 0. The fitted parameter 
values are listed in Table 2. Consistent with the results 
of the previous section, |Fi|/|i^2| ~ 2 at TOOK, and is 
almost independent of temperature down to 290 K. Sim- 
ilarly, rii, 0,2, and A are also independent of tempera- 
ture. On the other hand, F2 (the TA mode linewidth) 
is larger at 500 K and 290 K, which is consistent with 
the TA linewidth broadening below Td reported in Ref. T 
The linewidth of the TO mode Fi is also slightly broader 
below TOOK. 

The right panel of Fig. 5 shows the temperature de- 
pendence of the TA phonon for Q = (3,0,-0.2) mea- 
sured at (d) 5T0K, (e) 350 K, and (f) 50 K. Because 
these data do not include the TO mode, the profiles 
for 5T0K and 350K were fit fixing Qi ^ 11.2, Fi = 6, 
and A = IT. 4, while F1F2 > 0. Consequently, the un- 
certainties associated with the parameters in this zone 
will be larger. Values for the fitted parameters are listed 
in Table 2. From this analysis we find that the ratio 
|Fi(200)|/|Fi(300)| w 1. Since the nuclear structure fac- 
tor calculations imply |i^2(200)|/|i^2(300)| « 10, and we 
know that |Fi(200)|/|i^2(200)| « 2, we can estimate that 
|i^i(300)|/|F2(300)| « 20. This agrees reasonably well 
with the fitted ratio |Fi|/|F2| ~ 16 at both tempera- 
tures, and represents a reassuring self consistency check 
of our model. Moreover, a broadening of the TA mode is 
found at (3,0,-0.2) that is qualitatively similar to that 
observed at (2,0,-0.2). Thus, the temperature depen- 
dence of the phonon profiles is consistent with the mode- 
coupling model. However we did not obtain a clear TA 
phonon peak at 50 K. Instead, the solid line shown in 
Fig. 5(f) represents the profile that would be expected 
for the TA mode assuming the same parameters as those 
obtained from fits to the 5T0K data. The agreement be- 
tween the data and fit in this case is quite good. 
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FIG. 5: The left panel shows phonon profiles measured at Q 
= (2, 0, -0.2) at (a) TOOK, (b) 500 K, and (c) 290 K. The right 
panel shows TA phonon profiles measured at Q = (3, 0, —0.2) 
at (d) 570 K, (e) 350 K, and (f) 50 K. The Unes, except those 
shown in (f), represent fits to the resolution-convoluted MC 
function plus a Gaussian function centered at Tiuj = 0. The 
solid and dashed lines correspond to the MC and Gaussian 
cross sections, respectively. The solid line in (f) is an MC 
cross section using the same parameters as those obtained in 
(d). 



C. Dynamic structure factor for TO cross section 

The ratio of the dynamic structure factors that govern 
the transverse optic phonon intensities in different Bril- 
louin zones depends on the associated atomic vibrational 
displacements. Two of the most important atomic vi- 
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brational modes are those proposedJay Slater and Last, 
and are discussed by Harada et al£3 The Slater mode 
corresponds to atomic motions in which the oxygen and 
Mg/Nb atoms move in opposition while the Pb atoms re- 
main stationary. The Last mode corresponds to opposing 
motions of the Pb atoms and rigid (Mg/Nb) Oe octahe- 
dra. The dynamic structure factor is given by 

F^nei^Y.^Q-^,he^''-'''e~'^', (6) 

3 

where is the normalized displacement vector for the 
jth atom. Assuming only the Slater and Last modes are 
significant, then can be written in terms of the atomic 
displacement vectors for the Slater (sji) and Last (sj2) 
modes, and the parameter S, which is a measure of the 
relative contributions of the Last and Slater modes, as 
follows: 

= Sji + Ssj2. (7) 

For PMN, Hirota et aZ.0 calculated the value of \Finei\'^ 
for -0.5 < S* < 2 at (200), (300), and (110) as shown in 
Fig. 6 of Rcf. 11. Our result of |Fi(200)|/|Fi(300)| « 1 
obtained from our fits to the MC function is consistent 
with the value S' = 1.5 determined by Hirota et al. As 
an additional check, we measured the TO cross section 
at 80 K at several zone centers in the [hhl] zone. Be- 
cause the TO cross section at 80 K is greatly dimin- 
ished by the Bose factor, we were unable to perform 
a quantitative comparison of the phonon intensities to 
the dynamic structure factors using Eq. ^. However, we 
find that the TO phonon cross section at (111) is quite 
small. This observation agpes qualitatively with the cal- 
culations of Hirota et aZ.Ej using S = 1.5 for which 
|F„e^(lll)lVl^^ne^(200)|2 = 0.14. 



V. DISCUSSION 

In the previous section we presented our mode- 
coupling analysis and demonstrated that the character- 
istic features of the phonon profiles in PMN can be ex- 
plained by a coupling between the TO and TA modes, 
without need for any additional mode. This fact im- 
plies that the TO mode we observe is the ferroelectric 
soft mode. However our interpretation canfiicts with the 
model proposed by Naberezhnov et aL,H and later re- 
fined by Vakhrushev and ShapiroJiH which requires an 
additional soft mode. To reconcile these two models, we 
made a careful comparison of their data with ours. We 
find that, where there is overlap, the data agree quite 
well, and that only the interpretations are different. In 
this section, we attempt to explain the possible reasons 
for this difference. 



A. Naberezhnov- Vakhrushev branch 

In 1999 Naberezhnov et al. published an extensive 
study of the lattice dynamics of PMNE They observed 
well-defined TA and low-lying TO modes in the vicin- 
ity of the (220) zone center at 800 K. Surprisingly, they 
found that the structure factors of the TO modes are 
inconsistent with those derived from the diffuse scattei:^. 
ing intensities reported earlier by Vakhrushev et a/B'E£l 
This conclusion was based on the fact that the normal- 
ized displacements for the Pb {5{Ph) = 1.00), Mg or Nb 
{5{MN) = 0.18), and oxygen atoms {5{0) = -0.64), 
determined by Vakhrushev et aZ.,E3 do not satisfy the 
center-of-mass condition. Naberezhnov et al. noted that 
the center of mass of the unit cell is shifted from its 
original position, and thus cannot correspond to optical 
lattice vibrations. They attributed this shift to the slow 
relaxation of superparaelectric clusters. They therefore 
concluded that the observed low-lying TO mode could 
not be the ferroelectric soft mode. Instead they identified 
it as a hard TOl mode. A "quasi-optic" (QO) was mode 
was derived by fitting the TA profiles at the (200) and 
(110) zones to a function that couples the TA and QO 
modes, and this QO mode was identified as the soft mode. 
However, the inconsistency between the structure factors 
associated with the diffuse scattering and the TO mode 
has since been resolved by the "phase-shifted condensed 
soft mode" model proposed by Hirota et a^.llil Hirota et al. 
realized that the atomic displacements can be separated 
into two components: a uniform phase shift 5 shift — 0.58, 
and displacements that do satisfy the center-of-mass con- 
dition, namely 5cm{Ph) = 0.42, 6cm{MN) = -0.40, and 
5cm{0) = —1.22. Moreover, these center-of-mass dis- 
placements are consistent with the observed phonon cross 
section in that they correspond to a TO mode that is a 
mixture of Slater and Last modes in the ratio S = 1.5. 
Thus, the "phase-shifted condensed soft mode" model is 
able to reconcile the intensities of the TO phonon and dif- 
fuse scattering, and demonstrates that the diffuse scat- 
tering originates from the condensation of the soft TO 
mode. 

Vakhrushev and Shapiro have since analyzed the 
phonon profiles in the (300) zone where the TA phonon 
cross section is expected to be extremely weak as a result 
of a small dynamic structure factor. Their analysis uses 
a function that couples the TA and QO modes and ne- 
glectSrthe TA cross section because it was expected to be 
weak.tS But as demonstrated in Fig. 3(b) and in the ap- 
pendix, the coupling between modes greatly enhances the 
TA phonon cross section in spite of the small structure 
factor. Thus the TA mode is ncjt,negligible. The data 
taken by Vakhrushev and Shapirolij shows peaks that are 
completely identical to the TA peaks we observe. For ex- 
ample, they observe a peak at « 2 meV at (3,0.075,0) at 
880 K (Fig. 3 in Ref. 17), which is completely consistent 
with the TA peak we observe in Fig. 4(b). Another data 
of the identical sample at (3,-0.1,0) at 880 K is also 
published in Fig. 3 of Ref. 4. However, these peaks were 
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FIG. 6: Temperature dependence of the inverse integrated 
intensity of the zone center TO phonon. The inset shows 
TO phonon profiles measured at 1100 K and 700 K. The solid 
and dashed lines correspond to a damped harmonic oscillator 
description of the TO mode and a Gaussian form centered 
on hujo = 0, respectively (see text). The inverse integrated 
intensities are calculated from a numerical summation of the 
TO spectral weight corrected for the Bose factor. The high 
temperature data (T > 800 K) agree with the Gurie- Weiss 
behavior measured by Viehland et al, shown by the dotted 
line, for which 1/e oc (T - To), where To = 400 K. 



treated as a "Bragg tail"E3 in the analyses of Ref. 17. 
Again, the TA cross section is enhanced by the coupling 
and therefore one cannot neglect the TA cross section 
even though the structure factor for the TA mode in the 
(300) Brillouin zone is very small. Our model of coupling 
between TA and TO modes describes all of the charac- 
teristics of the phonon profiles and thus removes the need 
for an additional QO branch. 



B. Dielectric constant 

Because the QO mode is an artifact of the mode- 
coupling cross section, the Curie- Weiss behavior of the 
static dielectric constant must come from the softening 
of the zone-center TO mode. Therefore this TO mode 
is the ferroelectric soft mode. In fact, the zone-center 



TO mode measured at (200) exhibits a significant tem- 
perature dependence as shown in Fig. 1, and the linear 
temperature dependence of {hivo)'^ shown in Fig. 1(b) 
is qualitatively similar to that observed jiii conventional 
soft-mode ferroelectrics, such as PbTiOs.EJ However, the 
temperature dependence of {hujo)'^ in Fig. 1(b) does not 
satisfy the Lyddane-Sachs- Teller (LST) relation 



1/e oc {hujo)'^ 



(8) 



where e represents the static dielectric constant. Mea- 
surements of e by Viehland et alrn for T > 600 K show 
that the static dielectric constant exhibits a Curie- Weiss 
behavior 



ecxl/(r-To)^ 



(9) 



with To = 400 K, which is shown by the dotted line in 
Fig. 1(b). On the other hand, a linear extrapolation of 
{fiLUoY for T > Td gives Tq « OK. This was pointed out 
by Vakhrushev and Shapiro in Ref. 17. We have since 
reanalyzed the data and confirmed that the temperature 
dependence of {Tiojq)'^ indeed follows the solid line shown 
in Fig. 1(b). 

Vakhrushev and Shapiro found that T//(T) for T > Ta 
varies as T - Tq with Tq = 340 K, where I{T) is the in- 
tegrated intensity measured at (3, 0, 0) over the energy 
range 2 < ?iuj < 12meV. We have performed a similar 
analysis on our data measured at the (200) zone center 
which shows well-defined TO peaks. We fit the (200) zone 
center phonon profiles to a resolution-convoluted cross 
section composed of a Gaussian component to describe 
the elastic {hcu = 0) scattering plus a damped harmonic 
oscillator function to describe the TO mode. The fits to 
the data at HOOK and 700 K are shown in the insets to 
Fig. 6. The solid lines represent the damped harmonic 
oscillator component. We next evaluated the Bose- factor 
corrected integrated intensity for the damped harmonic 
oscillator component by numerical summation. Figure 6 
shows the temperature dependence of the inverse inte- 
grated intensity so obtained. The high temperature data 
agree with the Curie- Weiss behavior of 1/e reported by 
Viehland et al., with Tq = 400 K, which is shown as a 
dotted line in Fig. 6. 



C. Concluding remarks 

We have established the coupled nature of the TA and 
TO modes in PMN at 690 K, which lies above the Burns 
temperature Td- This coupling, along with the soft char- 
acter of the TO mode, leads naturally to the idea of a 
soft coupled optic mode that condenses at Td, and which 
therefore contains a significant transverse acoustic com- 
poncnL.^ This idea was originally suggested by Y. Ya- 
mada.o The condensation of such a coupled optic mode 
would then provide an elegant explanation for the origin 
of the phase-shifted condensed-soft mode model of the 
PNR proposed by Hirota et alM If correct, this concept 
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FIG. 7: Mode-coupling profiles as a function of Fi with (a) 
FxF2 = -2 (negative), and (b) F1F2 = 85 (positive). The 
other parameters were held fixed at the values Hi = 11, f22 = 
4, F2 = 1, and A = 20, which correspond to those obtained 
for (2,0, -0.2) and (3,0, -0.2). For negative Fifi, the spec- 
tral weight of the damped TO mode shifts toward the region 
between the TA and TO modes, and causes the apparent TA 
peak position to move to higher frequency. On the other hand, 
for positive F1F2, the spectral weight tends to shift towards 
the elastic (hu) — 0) position, and this causes the apparent 
TA peak position to move to lower frequency. 



may be the key to understanding the fundamental mech- 
anism underlying the unusual relaxor behavior. We have 
not as yet established a definite connection between such 
a soft coupled TO mode and the phase-shifted condensed 
soft mode. However this is one of our future projects. 
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FIG. 8: Mode-coupling profiles as a function of A. Values for 
the other parameters were held fixed at fii = 10, ^2 = 3, Fi = 
8, F2 = 1, and F1/F2 — 24, which correspond approximately 
to those obtained at Q = (3,0,-0.12) in Fig. 4(a). Part 
of the TO phonon intensity is shown to transfer to the TA 
mode as a result of the coupling, and results in a substantial 
enhancement of the TA mode intensity. 



VI. APPENDIX 

The phonon cross sections presented here have all been 
fit using the MC function (Eqs. 2, 3, and 4) given by 
Harada et al23 In this appendix we present several im- 
portant features of the MC function that give rise to large 
differences in the phonon profiles measured in different 
Brillouin zones. 

The effects of coupling can vary substantially when 
the sign of F1F2 is changed. Figure 7 shows the cross 
sections calculated using the MC function. The prod- 
uct of the two dynamic structure factors F1F2 was set 
to be —2 (negative) for Fig. 7(a), and 85 (positive) for 
Fig. 7(b). In both figures three cross sections are plot- 
ted for Fi = 12, 8 and 4, while the other parameters 
are fixed at T = 690 K, fli = llmeV, fl2 = 4meV, 
F2 = 1 meV, and A = 20. These parameters are approx- 
imately the same as those for obtained from fits to data 
measured for Q = (2, 0, -0.2) and (3, 0, -0.2). The inten- 
sities have been normalized so that the integrated cross 
section is kept constant. The TO phonon cross section 
for negative F1F2 shifts increasingly to the energy region 
between the TA and TO modes with increasing Fi. On 
the other hand, the TO phonon spectral weight shifts to 
the energy region below the TA mode near huj = for 
positive FiF2. Because of this, the apparent TA peak po- 
sition tends to shift to higher energy for negative F1F2, 
and to lower energy for positive FiF2. Thus, a change in 
the sign of F1F2 causes the apparent TA peak position 
to shift in different directions. 

Another important effect of the mode-coupling is 
shown in Fig. 8, which displays MC cross sections for 
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A = 0, 5, and 10. The other parameters are fixed at 
T = 690K, = 10, = 3, Ti = 8, T2 = 1, and 
F1/F2 — 24. Again the cross sections have been normal- 
ized to keep the integrated cross section constant. As 



a result of the coupling between modes, the TO cross 
section is partially transferred to the TA cross section. 
With sufficiently large values of A, a remarkably large 
enhancement of the TA peak can occur. 
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